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Abstract 

An extensive study of three-nucleon force effects in the entire phase space of 
the nucleon-deuteron breakup process, for energies from above the deuteron 
breakup threshold up to 200 MeV, has been performed. 3N Faddeev equations 
have been solved rigorously using the modern high precision nucleon-nucleon 
potentials AV18, CD Bonn, Nijm I, II and Nijm 93, and also adding 3N forces. 
We compare predictions for cross sections and various polarization observables 
when NN forces are used alone or when the 27r-exchange Tucson-Melbourne 
3NF was combined with each of them. In addition AV18 was combined with 
the Urbana IX 3NF and CD Bonn with the TM' 3NF, which is a modified 
version of the TM 3NF, more consistent with chiral symmetry. Large but 
generally model dependent 3NF effects have been found in certain breakup 
configurations, especially at the higher energies, both for cross sections and 
spin observables. These results demonstrate the usefulness of the kinemati- 
cally complete breakup reaction in testing the proper structure of 3N forces. 

21.30.-x, 21.45,+v, 25.10.+S, 24.70.+S 
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I. INTRODUCTION 



The nucleon induced deuteron breakup reaction together with elastic nucleon-deuteron 
(Nd) scattering have been considered for long as a valuable tool to test the three-nucleon 
(3N) Hamiltonian, in particular to shed light on the importance and structure of a three 
nucleon force (3NF) in the potential energy of the 3N system fl]]^]. However, up to recently 
difficulties in solving rigorously 3N equations in the continuum for realistic forces prohibited 
clear statements in such a study. The rapid progress in supercomputer technology made it 
now possible to gain numerically exact solutions of 3N Faddeev equations for any nucleon- 
nucleon (NN) interaction, even including a 3NF [Q. 

Present day models of NN interactions, mostly phenomenological and/or based on a me- 
son exchange picture, have achieved a high degree of maturity. It was possible by adjusting 
their parameters to describe with high precision (x 2 /data ~ 1) NN data in the large energy 
range from threshold to about 350 MeV |Q-|7J . These so called realistic potentials AV18 j5| , 
CD Bonn |||7|], Nijm I, II, and Nijm 93 [|J, are now extensively used in few-nucleon studies. 
Recent developments of chiral perturbation theory lead to NN potentials, where the one- 
and multi-pion exchanges are treated unambiguously according to chiral symmetry 



Therefore the theoretical uncertainties of meson exchange based models have been at least 
partially removed and this step puts the theory of nuclear forces on a more solid basis. First 
successful applications of these interactions in 3N and 4N systems have already demon- 
strated their power in interpreting and describing the data [|Hj-|2~T||. In this paper, however, 
we shall only apply the above mentioned high-precision NN interactions. 

This present state of NN potential models, describing NN data perfectly, suggests now, 
when going to systems with more than two nucleons, to concentrate on the significance and 
properties of many-body force contributions to the nuclear Hamiltonian. In case of three 
nucleons the first indication on the importance of 3NF contributions came from 3 He and 3 H 
bound state studies. All realistic NN potentials are unable to reproduce their experimental 
binding energies leading to a 3 He and 3 H underbinding of the order of 0.5 -j- 1 MeV |[22|| . 
This clear cut underbinding exists also for 4 He where it amounts to 2 -=-4 MeV |22-23]. Also 
for higher mass nuclei up to A=8, where stochastic techniques must be applied, realistic NN 
forces failed to provide the experimental binding energies fl2~5|j26|| . 

A natural step to explain this underbinding was to consider 3NFs in the nuclear Hamil- 
tonian. Presently, the most often used dynamical process is the 27r-exchange between three 
nucleons. An important dynamical ingredient in that process is, as suggested long time ago 
by the Fujita-Miyazawa 3NF [p7| , an intermediate excited nucleon state A. It was augmented 
later by further ingredients leading to the Urbana IX p8| and Tucson-Melbourne (TM) |29| 



3NFs, which are mostly used in present day few-nucleon calculations. It was possible by 
properly adjusting the parameters of these 3NFs to get essentially correct 3N and 4N bind- 
ing energies. Also a fairly successful description of low energy bound state energies of up 
to A=8 nuclei resulted when adding 3NFs. This was recently improved by adding further 
3NFs related to three-pion exchange with intermediate A's |R0 . In the present exploratory 
investigation we shall not yet include this newest generation of 3NFs. 

Though this first signal on 3NF effects, resulting from discrete states, is important an 
even more sensitive and detailed investigation of 3NF properties can be carried through 
in scattering processes, where already a rich set of spin observables has been measured 
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and further data are expected to come pT|-|3~4}j . From the theoretical side one can make 
exact predictions for these observables using nuclear forces in all their complexities ||. 
Experimentally one can access nowadays spin observables in elastic Nd scattering where in 
the initial states the deuteron and/or the nucleon is polarized and where in the final state 
the polarizations of the outgoing particles can also be measured |35|-f43"|| . This together with 
cross sections leads to a very rich spectrum of observables in Nd elastic scattering and the Nd 
breakup process. Such a set of observables will be a solid basis to test the 3N Hamiltonian. 
Using available model Hamiltonians one can provide guidance to select specific observables 
and energies which are most appropriate to see 3NF properties. In [33 such a study has 
been performed for elastic Nd scattering. It was shown that the current 3NF models exhibit 
a lot of effects and more data are needed to provide constraints on the theoretical models 
of 3NFs. First data sets pl| , |35|j38| -^0| in elastic Nd scattering showed both, successes and 
failures of present day 3NF models combined with the realistic NN forces. It is the aim 
of this paper to probe the cross section and several spin observables of the kinematically 
complete Nd breakup process against sensitivity to 3NFs combined with the high-precision 
NN interactions. In a following paper (II) already existing data will be compared to those 
model predictions. 

For the convenience of the reader we briefly review in Section [TT] our theoretical formalism 
for the 3N continuum and display the 3NF models which are presently in vogue and which 
we use. Our predictions for quite a few breakup observables based on various nuclear force 
combinations are presented in Section |ITJ While in this Section we investigate the full 
3N breakup phase space we concentrate on a subset of breakup configurations in Section [TV 



according to the experimental restrictions imposed in an ongoing experiment at KVI p3| , f45 
We conclude in Section [V[ 



II. THEORETICAL FORMALISM AND 3NF MODELS 

We refer to H for a general overview on 3N scattering and specifically for our way to 
formulate it. The inclusion of 3NFs follows [f|(J. This is a direct generalization of what is 



being used for the 3N bound state JTj]]. We define an amplitude T via our central Faddeev- 
like equation: 

T = tP4> + (1 + tC ) U 4 (1) (l + P)4> + tPG T + (l + tG )v} 1) (1 + P)G T (1) 

The initial channel state is composed of a deuteron and a momentum eigenstate of the 
projectile nucleon. The two-body t-operator is denoted by t, the free 3N propagator by Go 
and P is the sum of a cyclical and an anticyclical permutation of three particles. The 3N 
force V4 can always be decomposed into a sum of three parts: 

V, = v^ + Vi 2) + Vf\ (2) 

where V® is symmetrical under the exchange of the nucleons jk with j 7^ % 7^ k. As seen 
in Eq. (|l|) only one of the three parts occurs explicitly, the others via the permutations 
contained in P. 

The physical breakup amplitude is given via 
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U = (1 + P)T. 



(3) 



The Faddeev-like integral equation ([!]) has the nice property that the terms resulting by 
iteration and inserted into Eq. (|D yield the multiple scattering series, which gives transparent 
insight into the reaction mechanism. 



The definition of the various spin observables can be found in p j4q , |49|l . We shall en- 
counter nucleon and deuteron vector analysing powers A y and A^\ where in the initial 
state either the nucleon or the deuteron is polarized. Further, the deuteron can be tensor 
polarized in the initial state leading to the three tensor analysing powers T 2 k (k = 0, 1, 2) or 
the corresponding cartesian analyzing powers Ay = x,y, z). 

We concentrate in this study on the kinematically complete 3N breakup. The final 3N 
state of this reaction is specified by nine momentum components. In such an experiment 
the energies and directions (polar and azimuthal angles 9i and fa) of two outgoing nucleons 
(i = 1, 2) are measured. This together with restrictions imposed by energy and momentum 
conservation lead to a correlation of the measured energies E\ and E 2 , causing them to lie 
on a kinematical curve (S-curve) in the Ei — E 2 plane. All observables will be presented as 
a function of the arc length S related to a point along this kinematical curve. It is a matter 
of convention to choose the location of S= on the S-curve and we use the one defined in 
Ref. f. 

As NN forces we use the five realistic NN interactions mentioned in the introduction. 



They are combined with 3NF models. The 27r-exchange TM model |29] has been around for 
quite some time and is based on a low momentum expansion of the 7r — N off- (the mass-) 
shell scattering amplitude. It has the following form: 

Vi 1] = % H(Q 2 )H (Q' 2 



x jf a • r 3 (a + bQ ■ Q' + c(Q 2 + Q' 2 )^j + d ir 3 x r 2 ■ f x B x ■ Q x Q'j . (4) 

The elements of the underlying Feynman diagram are obvious: the two-pion propagators 
depending on the pion momenta Q and Q', the strong form factors combined into the two H- 
functions and most importantly the parameterization of the irN amplitude inside the curly 
bracket which is combined with the isospins T2 and T3 of the two accompanying nucleons. 
The if -functions are parameterized as: 



A 2 - ml 

,A 2 + g s 



(5) 



In what we denote by the TM 3NF we use the original parameters a = 1.13/77^, b = 
-2.58/m;[, c = 1.0/ml, d = -0.753/m^ where m n = 130.6 MeV @. The still very 
demanding computer resources prevented us to use the recently updated values |51|] which, 
however, differ only slightly from the one we use and we do not expect significant changes. 
The b- and c/-terms are mostly determined by an intermediate A in a static approximation. 
The a- and c-terms are related to 5-wave 7rN scattering. The cut-off parameter A is used 
to adjust the 3 H binding energy separately for different NN forces [22]|. In units of the pion 
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mass we find the A to be: 4.856, 5.215, 5.120, 5.072, and 5.212 when the TM 3NF is 
combined with CD Bonn, AV18, Nijm I, II, and Nijm 93, respectively. 

Of course in a meson exchange picture additional processes should be added containing 
different mesons like n — p, p — p, etc.; different intermediate excited states might also 
play a role |5"2fl . To some extent 3NF models with respect to those extensions have already 
been developed and applied [^3] -|57|]. Our adjustment of A is a very rough manner to take 



other processes into account. Further studies incorporating these additional 3NFs should be 
performed both for elastic Nd scattering and the breakup process. 

The parameterization of the TM 3NF has been criticized, since it violates chiral symme- 
try [p^y^l . A form consistent with chiral symmetry (though still incomplete to that order 
in the appropriate power counting) is obtained by keeping the long range part of the c-term 
leading to a changed parameter a' = a — 2m 2 = — 0.87/77^ [|5£|[5£j] and dropping the short 
range part. This essentially means a change of sign for a. This form will be called TM' 
later on. The corresponding A-value, when TM' is used with the CD Bonn potential, is 
A = 4.593. 

The two-meson exchange 3NF has been also studied by Robilotta et al. 
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to the Brazilean 3NF. It is similar to the one of TM and also the results gained for 3N 
observables |6T|] are similar to the ones for the TM 3NF. In this paper we do not take 
that force into account. Instead we included the Urbana IX 3NF [PSfl, which is intensively 



used in the Urbana- Argonne collaboration. That force is based on the old Fujita-Miyazawa 



ansatz [27[ of an intermediate A occurring in the two-pion exchange and is augmented by a 
spin and isospin independent short range piece. It has the form 



(i) 



.4 



2tt 



{Xi2, X 13 } {ft • r 2 , n ■ f 3 } + i [Xi2, X 13 ] [ft ■ f 2 , fi ■ r 3 



+ 



+ U T*(r 12 )T*(r 13 ), (6) 

where {, } and [, ] are anticommutator and commutator brackets and X^ is defined as: 

Xij = Y^rij) Bi ■ ffj + T n (rij) Sij (7) 

Here Oi is a nucleon spin operator, is the standard tensor force, the T(r i:) ) and Y(rij) are 
functions of the distance nj between nucleons i and j, that are associated with the Yukawa 
and tensor part of the one-pion-exchange interaction 



YJr) 



1 - e" 



(8) 



1 + 



(9) 



The parameters 



quenched via a short-range cut-off function (1 — e _cr ) with c=2.1 fm~ 2 |62 
A 2n and Uq were chosen to reproduce the experimental 3 H binding energy and to obtain the 
empirical equilibrium density of nuclear matter when this 3NF model is used with the AV18 
potential. Their values for model IX are A 2w = -0.0293 MeV and U = 0.0048 MeV J28 
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Since we work in momentum space using a partial wave expansion, the form given in Eq. 



has to be rewritten. This has been done in 44 



Since there is no apparent consistency of the mostly phenomenological realistic NN forces 
and the 3NF models, we test various combinations thereof. In all cases, however, we require 
that the particular force combination should reproduce the experimental triton binding 
energy. Some of the 3N observables scale with the triton binding energy |63|]. The adjustment 
to the triton binding energy has the advantage that our investigation is not misled by these 
scaling effects. 

With respect to the intermediate A one should say that very likely the static approxi- 
mation is not justified and the A should be allowed to propagate like the nucleon. This has 
been pursued intensively, for instance, by the Hannover group [64]] and their recent work has 
been also devoted to the 3N continuum [ |65| . 

In view of all that, it is quite clear that our present study is not at all complete but can at 
least provide some insight into effects specific 3NF models generate for breakup observables. 



III. PREDICTIONS OF 3NF EFFECTS 

Since we would like to cover a wide range of incoming nucleon energies, ranging from just 
above the Nd breakup threshold up to 200 MeV, it is necessary to take a sufficient number of 
partial waves into account to guarantee converged solutions of the Faddeev equations. In all 
presented calculations we went up to the two-nucleon subsystem total angular momentum 
jmax = 5. This corresponds to a maximal number of 142 partial wave states (channels) in 
the 3N system. We checked that the convergence has been achieved by looking at the results 
obtained for j max = 6, which increases the number of channels to 194. This convergence 
check refers to a calculation without a 3NF. The inclusion of 3NFs has been carried through 
for all total angular momenta of the 3N system up to J = 13/2. These high angular momenta 
are required at the higher energies > 100 MeV. The longer ranged 2N interactions require 
states up to J = 25/2 at the higher energies in order to get converged results. 

A phenomenological criterion for 3NF effects is that they lie outside of the spread of 
the five realistic NN force predictions only. This spread will be indicated by a band (called 
the "2N" band). Unfortunately we cannot include the pp Coulomb force effects. At the 
higher energies, however, those effects should be small. Also in case of AV18 we do not 
take the various electromagnetic corrections into account, which leads among possibly other 
effects to a slightly wrong deuteron binding energy (Eb = 2.242 instead of 2.225 MeV). 
This shifted ^-value, however, has only a small effect on our results, which is mostly of 
kinematical origin. Also the NN-phase shifts, obtained without those additional terms, differ 
only slightly from the standard ones. The kinematical effect also leads to a slight shift in 
the position of the S-curve. 

To investigate 3NF effects we combine the TM 3NF with the five NN forces forming 
a second band (called the "2N+TM" band). In all cases the cut-off value A in Eq. @ 
has been adjusted separately for each NN force to the 3 H binding energy |66| as given in 
the previous Section. Since that interplay is a purely phenomenological step, the outcome 
is theoretically not under control and we combine all the results into a second band. In 
addition, we want to compare the TM 3NF and the modified TM', which is more consistent 
with chiral symmetry. We combine TM' with CD Bonn and display the results as a dashed 
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curve. Finally we compare the AV18 with the Urbana IX 3NF what is presented as a solid 
line. There are clearly more combinations possible but we felt this to be sufficient to get an 
orientation on the magnitudes of the effects. 

Before starting to investigate exclusive breakup configurations with large 3NF effects, 
it is interesting to look at the elastic scattering and breakup contributions to the total 
cross section for neutron-deuteron (nd) scattering. Fig. [I] shows the total cross section 
atot, predicted by the CD Bonn potential, as a function of the incoming neutron laboratory 
energy up to 300 MeV together with the contributions coming from elastic scattering and 
the breakup processes. It is seen that the total elastic cross section decreases rapidly with 
increasing energy of the incoming neutron. At about 50 MeV the breakup contributes 
already ~50% to the total cross section. At higher energies it dominates the total cross 
section and at about 200 MeV the elastic scattering contributes only ~20% to the total 
cross section. This picture does not depend on the choice of the NN potential used: the 
theoretical results for the nd total elastic and breakup cross sections are stable under the 
exchange of the realistic NN potentials. The differences are of the order of 1%. 

In Figs. we compare the theoretical predictions with the experimental data. There is 
a good agreement between the "pure" NN theory and the data up to ~50 MeV for all three 
total cross sections. At higher energies the data exist only for a to t and the theory begins to 



underestimate the data (about 2% at around 100 MeV and about 13% at ~300 MeV) |67] , |68 

In order to see the influence of 3NFs, we calculated the total cross section for nd scattering 
combining some of the NN potentials and 3NF models as described above. We see, adding 
a 3NF increases the value of the different total cross sections and shifts the theoretical 
predictions closer to the data (it is clearly visible for Ei a b > 100 MeV in Fig. |5|). The 
predictions obtained with different 3NF models are practically the same. 

In order to quantify the effect of a 3NF on the total cross sections we calculated its 
relative contributions in the corresponding processes as: 



3NF 



<2N) (2N+TM) 

07 -a) 



x 100% (10) 



where <Tj = a tat , uf ^ s , cr b t l t . As can be seen from Fig. the 3NFs enhance the nd total cross 
section only by about up to 5%. This enhancement is practically independent from the 
particular 3NF used and for energies above 100 MeV changes only slightly with energy. Due 
to the drastically diminishing contribution of the total elastic cross section most of the 3NF 
effects at higher energies must occur in the breakup process. Since the total breakup cross 
section is dominated by the quasi free scattering, for which 3NF effects are negligible P,|69|, 
it follows, that at higher energies 3NF effects must be found in other breakup configurations, 
possibly just in some specific geometries. It is one of the aims of this paper to investigate 
this conclusion in detail. 

As can be seen from Fig. [5], addition of a 3NF only partially removes the discrepancy 
between data and theory. At higher energies one has to expect that relativistic effects might 
appear. A generally accepted framework for carrying out relativistic 3N calculations does 
not yet exist. An estimation of these effects due to proper relativistic treatment of the 
kinematics while keeping the (non-relativistic) scattering amplitudes unchanged was done 
in fl68f . It was found, that such relativistic effects enhance the nd total cross section by 
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about 7% at 250 MeV while being much smaller at lower energies. These corrections bring 
the theoretical predictions for the nd total cross section much closer to the data. 

In the following we restrict the study of 3NF effects to the nd breakup reaction. Contrary 
to elastic scattering, where the description of the two-body final state requires only one 
variable (e.g. the cm. angle 9), the 3N final state of the breakup process is kinematically 
much more complex and requires five parameters to be described completely. The possibility 
to choose very many specific kinematical geometries and energy distributions for the three 
outgoing nucleons allows to probe the 3N Hamiltonian in greater detail. 

Our goal is to scan the full breakup phase space for the configurations which exhibit large 
3NF effects in some observables. In order to cover the broad range of incoming energies, we 
took, as representatives, the following laboratory energies: 13, 65, 135, and 200 MeV. 

As a first exploratory study, we calculated the five-fold differential cross section a = 
dsdfiidn 2 ' and tlie ana ly zm S powers A y , A x , A z , A y d \ A yy , A xz , A xx using the CD Bonn NN 
potential with and without the TM 3NF. The observables were tabulated as a function of 
Ei, E 2 , 0i, 9 2 and <pi 2 = <f>i — (f) 2 , varying 9{ in steps of 5°, 0i2 in steps of 10° and taking a 
step along the S-curve, A S(E 1 , E 2 ) = 1 MeV, over the full phase space. 

In order to check if the chosen grid over the full phase space was sufficiently dense we 
calculated: 

where N is the number of the calculated configurations. In Tab. | the values of cr^ n are 
compared with the total breakup cross section a\l t . The difference between a\l t and is 
smaller than 5%, what gives a strong support that our grid is sufficiently dense. 

In order to discover interesting cases we looked for configurations with the largest differ- 
ences between observables predicted with the CD Bonn alone and when the TM 3NF has 
been included. The influence of the TM 3NF can be quantitatively defined by the relative 
difference SO of the observable O = a, A y , A x , A z , A y d \ A yy , A xz , A xx at a given phase 
space point (Ei, E 2 , 9i, 9 2 , 4>i 2 ) when the TM 3NF is switched on 



SO(E 1 ,E 2 ,9 l ,9 2 , 



'12 j 



q( 2N ) - Q{2N+TM) 



Q(2N) 



x 100% (12) 



Here 0^ 2N ^ and (9( 2Ar +™) are the predictions for the particular observable O using the 
CD Bonn potential alone and including in addition the TM 3NF. In the case of analyzing 
powers, due to experimental limitations only kinematical points have been considered for 
which \Af N) \ > 0.05. 

In order to get insight into the limitations imposed by the experimental conditions in 
detecting interesting cases we performed the search in three steps. First, for each observable 
O at each energy we searched over the entire grid of the phase space for those points where 
50 achieves its maximum value AO: 

AO = max[60{E u E 2 , 6 U 9 2 , (f> 12 )} 

The resulting AO's are presented in the rows labeled all in Tab. |I[ Second, bearing in 
mind that we look for effects that should be experimentally accessible, we searched for the 
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maximum of 50 only at phase space points for which the cross section > 0.01 ?t rr 

and this for at least 5 MeV along the S-curve (a-cut). The resulting AO's are shown in the 
rows labeled er-cut. And finally, in addition to the er-cut, we demand that the energies of 
the detected nucleons (E\, E 2 ) must be greater than 15 MeVQ The final AO's are displayed 
in the third rows labeled a + £7-cut. Supplementary, in the last column of Tab. [II] labeled 
N we present the number of configurations fulfilling without maximum requirements the 
three search restrictions: no restrictions, applying the er-cut, and applying both, the a- and 
-E-cuts. 

The cross section (er-cut) and energy (E'-cut) threshold values are based on current 
experimental conditions and limitations. We are aware that each of these values is somewhat 
arbitrary. We do not need to lower them, however, since as we show below the 3NF effects 
are already clearly visible under these restricted conditions. 

As we can see in Tab. [Tl] the 3NF effects, in general, increase with the energy of the 
incoming nucleon. This tendency also holds after applying the er-cut and the er-cut together 
with the E-cut. The only exceptions are A x (where the maximum value of AA X including 
the cuts is at Ei ab = 65 MeV) and A y and A xx (where the 3NF effects are largest at 135 MeV 
rather than at 200 MeV). The values of AO after applying both the a- and E-cuts do not 
change in a considerable way but the number of configurations N decreases dramatically. 

The polarization observables are particularly sensitive to the 3NF. Among them the 
tensor analyzing powers are especially interesting, showing large 3NF effects. The nucleon 
and deuteron vector analyzing powers (A y , Ap) are good tools to test the 3NF as well. 
Predictions for the differential cross section are also very promising, especially at the energies 
135 and 200 MeV. The 3NF influence is the least evident in A z and A x . 

The next step in our analysis was to divide the phase space into regions according to 
certain values of 50 = SO(9\, 9 2 , <fii 2 , Ei, E 2 ) for our chosen observables. Of course, the 
regions where 50 takes its largest values are of the strongest interest. In order to display 
these regions, we project them onto three planes: 9\-6 2 , Q\-§\ 2 and E\-E 2 . Fig. [7] shows one 
example for Acr at 135 MeV. The three rows represent the three cut conditions described 
above and labeled as in Tab. [TJ. Using the information from the first and the second 
column for Acr, varying between and 90%, one can locate the angular regions in phase 
space where Acr values of a certain magnitude can be found. The last column provides the 
additional information in the energy distribution. For example our sample plot shows one 
very promising configuration ((6^ ~ 15°, 9 2 ~ 15°, (f>i 2 ~ 0°)) with Acr > 90%. Moreover that 
configuration survives all the applied cuts. The corresponding situation for A yy is shown in 
Fig. || It is seen that 3NF effects are larger in magnitude and they occur in different regions 
of the phase space compared to the interesting ones for Acr. 

In this way we found at each energy the most promising configurations for each studied 
observable. It turned out that the influence of the 3NF is not limited to separate points 
but rather shows up in extended regions of phase space. Also the regions differ with the 
observable chosen. 

Next we performed more detailed studies. For each observable and each energy we have 
chosen about 20 angular configurations and calculated the observable using the five high- 



x 3 MeV for E lab = 13 MeV 



9 



precision potentials. In the following figures they form the light shaded band. Then the TM 
3NF was combined with all five NN potentials forming a second band (dark shaded band in 
the following). In addition, we regard the force combinations CD Bonn with TM' and AV18 
with Urbana IX. 

As a result of our study, we have a huge amount of information on the 3NF effects 
for each observable at each of the four energies. In the following, we present only general 
conclusions about each observable and we only show the most interesting configurations. 
More results and examples can be found in |){| . 

A. 3NF effects at 13 MeV 

The case of 13 MeV serves as an example for 3NF effects in a low energy breakup process. 
From Tab. [IJ it follows that for vector analyzing powers no significant 3NF effects are present. 
The values of AO for tensor analyzing powers and the cross section give some hope to find 
configurations with measurable 3NF effects. In Fig. |9] we show some examples in which we 
found the largest effects for the differential cross section, A xx , A yy and A xz . Unfortunately, 
the spread in the predictions of various NN potentials can be as large as 10% and this is 
comparable in magnitude with the 3NF effects found. Therefore the conclusion is that the 
3NF effects at 13 MeV are too small to be distinguished from the dependences on the NN 
force model. 

B. 3NF effects at higher energies 

To look for the magnitude of 3NF effects at higher energies we investigated each of 
the another eight observables listed above at 65, 135 and 200 MeV. We focus only on 
configurations which survive both the a- and .E-cuts, because our aim is to find regions of 
the phase space with large and measurable 3NF effects. In all cases, however, we checked if 
the rejected configurations contain interesting 3NF effects. 

1. Differential cross section 

At higher energies there are significant effects of 3NFs in the differential cross section. 
They are particularly large (more than 90%) in some configurations at 135 and 200 MeV. 

Generally, in all phase space points studied adding the TM 3NF increases the value of 
the cross section at 135 and 200 MeV, whereas at 65 MeV we found configurations where 
the TM 3NF decreases a. We show prominent cases in Fig. 0. We see that at 135 and 200 
MeV 5a > 60% (in fact this happens not only for 12 = 0° but in a range as large as 20°). 
In the maxima of the cross sections the energies of two outgoing nucleons (Ei and E 2 ) are 
comparable and together amount to about 80% of the total accessible energy. The "2N" 
and "2N+TM" bands are well separated. 

Astonishingly, at 135 and 200 MeV the AV18+Urbana IX and CD Bonn TM' predictions, 
while close together, are distinctly different from all NN force combinations with TM. 
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2. Vector analyzing powers A y and Ay 

The 3NF effects for the nucleoli and deuteron vector analyzing powers A y and Af in- 
crease with energy and become quite large at higher energies in some specific configurations. 

Most of these configurations are close to the so-called final state interaction (FSI) geom- 
etry, where the two outgoing nucleons have equal momenta. For example, at 200 MeV the 
value of SA^f) reaches about 350% for the (9i, 8 2 , 0i 2 ) = (25°, 25°, 0°) configuration while in 
others i5ijf < 150%. Similarly, 5 Ay reaches the maximal value of nearly 200%, but only 
in few configurations goes above 140% (six for 135 MeV and three for 200 MeV). Figs. [Ll] 



and [L^ show examples of the most interesting configurations at each energy. Despite the 
substantial broadening of both "2N" and "2N+TM" bands (especially at 65 MeV), the TM 
3NF effects are clearly visible. At higher energies the predictions with and without TM 3NF 
have even different signs, both for A y and A y d \ We found also some configurations where 
the inclusion of the TM 3NF decreases the pure NN value of the analyzing powers from the 
positive value 0.2 to zero (e.g. (15°, 15°, 20°)). Data from such configurations can be an 
excellent tool to test 3NF models. 

We cannot draw any general conclusions about the predictions obtained using Urbana IX 
with AV18 and TM' with CD Bonn, since they are quite different in various configurations. 
In addition they can be also quite different from the TM predictions. 



3. Nucleon analyzing powers A x and A z 
To get nonzero values for the nucleon analyzing powers A x and A z one has to go out 



of the reaction plane |48]|49|]. For both of them the 3NF effects are the smallest among all 
considered analyzing powers. Moreover, the largest 5A X value was found at 65 MeV and it 
decreases with increasing energy (assuming cuts). For A z the maximum values of 5A Z at 
135 and 200 MeV are comparable (5A Z = 60%) and they are larger than the corresponding 
value at 65 MeV (5A Z =20%). For A z the assumed threshold value of \A Z \ (\A Z \ > 0.05) is 
crucial because \A Z \ reaches measurable values only at higher energies. There this cut is not 
so decisive anymore because the maximum of \A Z \ is larger and approaches 0.45. 



Again, we picked up the most promising configurations and show them in Figs. [L3| and 
14] for A z and A x , respectively. At all energies the "2N" and "2N+TM" bands are rather 
narrow and the 3NF effects are clearly visible. For A x and A z the Urbana IX and TM' forces 
provide smaller 3NF effects. 



4- Tensor analyzing powers 

These spin observables appear to be the most sensitive to the 3NFs we used. For each 
of them the magnitude of 3NF effects increases strongly with the energy. There are ex- 
tended regions of the phase space with large 3NF effects at each energy, what makes the 
investigation of A xx , A yy and A xz extremely promising. At the highest energy there are a 
few configurations with especially large values of 5Ai (Ai = A yy , A xx , A xz ). One distin- 
guished example is A yy for which SA yy reaches the extremely large value of about 400% 
in the (45°,45°,0°) configuration (see Fig. [Tj^). For this observable there are in addition 
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three configurations with 5A yy > 250%. They are close to FSI geometries with comparable 
energies of the two outgoing nucleons. We found a similar behavior for A xx and A xz . 

At the higher energies we even found several configurations with different signs for the 
"2N" and "2N+TM" bands what leads to a very big value of 8A4. However, they do not 
always represent the best region to investigate 3NF effects. For example, 5A XX at 200 MeV 
reaches the largest value for the (145°, 10°, 180°) configuration but the short length of the S- 
curve, where 3NF effects are large, makes this configuration hard to use for an experimental 
study. We have also found several configurations with large 3NF effects but the large spreads 
of "2N" and "2N+TM" bands inhibit the unambiguous conclusion about the significance of 
3NF effects, e.g. for A xz at 200 MeV in (45°,45°,0°) or configurations with largest 5A yy . 
On the other hand some of these configurations present themselves as excellent playground 
to look for differences between different models of NN interactions. 



Figs. [L5|, [16| and [L7| show our choice for the most sensitive configurations. The predictions 
using Urbana IX combined with AV18 and TM' with CD Bonn are often close to each other 
and give smaller effects than with the TM 3NF. 

IV. SEARCH FOR CONFIGURATIONS WITH LARGE 3NF EFFECTS IN THE 

P(D,PP)N BREAKUP AT 130 MEV 



Awaiting results from a new experiment currently under way at KVI in Groningen []33[ , we 
calculated the differential cross section dn ^n 2 ds together with vector and tensor analyzing 

powers A y d \ A xx , A yy , A xz for the p(d,pp)n breakup reaction at Ed = 130 MeV in the 
laboratory frame. All these observables are planned to be measured. The calculations were 
performed for 970 configurations covering the entire experimentally available phase space 
(the maximal available angle of the outgoing protons is 9 = 76.9° in the laboratory frame). 
We used the CD Bonn potential with and without the TM 3NF varying the angles 81, 62 
and <fti2 in steps of 5° and the position on the kinematical locus S(Ei, E2) in steps of 1 MeV. 
This exploratory study with the only force combination CD Bonn and TM was performed in 
the same way as in the previous Section |T| and the results are presented in Tab. [TTI[ Again 
the tensor analyzing powers are the most sensitive observables to 3NFs. Then we performed 
in addition a more detailed study using all force combinations as in the previous Section. 

A. Differential cross section 

In the case of the differential cross section the experimentally required condition Ei > 
15 MeV rejects a lot of very promising configurations with 12 < 70°. Despite that, we found 



many configurations with 3NF effects in the region of the maximal values (see Tab. 
Very interesting configurations are (61 ~ 20°, 62 — 20°, 0i2 > 130°) which correspond to the 
situation where one of the outgoing nucleons takes all accessible energy. 

In Fig. [18| we show cross sections for three configurations revealing one of the largest 3NF 
effects. Only part of the allowed S-curve is presented, for which the energies of the protons 
fulfill the conditions > 15 MeV. It is very encouraging that all of these configurations 
are accessible in the Groningen experiment. The NN force predictions are rather close 
together. Including the TM 3NF generally increases the cross section in these configurations. 
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Predictions based on the TM' and Urbana IX forces are of similar magnitude to the TM 
ones. 



B. Deuteron vector analyzing power 

The Atf) seems to be a good tool to investigate 3NF effects. The most sensitive config- 
urations occur where the energies of two outgoing protons are small (< 35 MeV). Unfor- 
tunately due to the experimental arrangement, which allows to measure only protons with 
Ei > 15 MeV many interesting configurations will not be studied in the experiment. Nev- 
ertheless we found extended regions of the phase space with 0i j2 > 20° and 0i 2 < 50° where 
the value of SA y d ^ reaches more than 60%. From these regions we picked a few configura- 
tions and show them in Fig. [L^. The "2N" and "2N+TM" bands are well separated showing 
clear and large (up to 70%) 3NF effects. The NN force band is rather narrow, whereas the 
"2N+TM" band is distinctly broader. For the chosen configurations the CD Bonn+TM' 
and AV18+Urbana IX are close to each other and predict very small 3NF effects, different 
from that of the TM 3NF. 



C. Tensor analyzing powers 

Finally, we consider the three tensor analyzing powers A xx , Ayy and A xz . The effects in 
A xx strongly depend on the energy and, unfortunately, the most interesting part of the phase 
space with largest effects is outside the range of the Groningen experiment (£1,2 < 15 MeV 
and 9i 7 2 > 40°). However, we found some regions where smaller but clear 3NF effects are 
visible (see Fig. ^). The predictions obtained using Urbana IX and TM' forces are close 
together and often they do not differ from the pure NN force predictions. 

For A yy we found that the most sensitive parts of the accessible phase space are for 
20° < 6\ < 25° and 9 2 < 25° with 12 < 60°. In Fig. [H| we show the three most interesting 
configurations. The "2N" and "2N+TM" bands are rather narrow and large 3NF effects 
rising above 100% are clearly visible. The Urbana IX and TM' predictions are rather close 
together and give smaller effects than the TM 3NF. 

For A xz the configurations with largest effects are for #i )2 = 20° and 0i 2 > 130°. These 
configurations correspond to the situation where the energy of one of the outgoing nucleons 



is much larger than the energy of the other nucleon. In Fig. |22] we show three examples 
where the 3NF effects are larger then 100%. Again the effects of Urbana IX and TM' are 
significantly smaller than the effects of TM 3NF. 

In conclusion, in the dN breakup reaction with an incoming deuteron energy of 130 MeV 
the 3NF effects are clearly visible for cross section, vector and tensor analyzing powers in 
numerous configurations. The most promising observables are the tensor analyzing powers 
where the 3NF effects are as large as 100-180% . For these observables the magnitudes of 
the effects depend on the 3NF model used. 
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V. SUMMARY AND CONCLUSIONS 



We started our investigation with the nd total cross section. As could be foreseen we 
found that the magnitude of 3NF effects increases with energy. Because at higher energies 
the breakup contribution to the total nd cross section becomes dominant this implies that 
the breakup process should be especially gratifying to see 3NF effects. In such a process, 
where the three final momenta and, in addition, also some spin orientations are fixed, one 
has the most sensitive insight into the 3N Hamiltonian. Therefore we performed a detailed 
study of the nd breakup reaction focusing our investigations on the effects of 3NFs on the 
cross section and quite a few polarization observables. We compared predictions based on 
realistic NN forces (CD Bonn, AV18, Nijm I, II and 93) to those obtained when, in addition, 
three different 3NF models (TM, TM' and Urbana IX) were included in the 3N Hamiltonian. 
The NN force predictions only have a certain spread. Also adding 3NFs to them yields again 
a spread. We talk of 3NF effects if the two spreads are clearly distinct. We have found in 
cross sections and analyzing powers large 3NF effects which clearly should be discernible by 
experiments. 

In our study we covered the entire phase space at several energies: 13, 65, 135 and 200 
MeV. We located for each observable the regions in phase space where the largest 3NF 
effects occur. These regions are in general different for the various observables. The 3NF 
effects increase in general with energy. Especially promising will be the tensor analyzing 
powers A xx , A yy and A xz . Also the five- fold differential cross section turned out to be a 
promising observable to clearly see 3NF effects. The smallest effects have been found for 
the analyzing powers A x and A z . We did not study all possible force combinations between 
the five high precision NN forces and the three 3NF models, which are currently most often 
used. Our study aimed only at a first orientation. But the results show that the outcome 
for the 3NF effects can depend very much on the force combination. This tells that the 
3NF effects predicted should be considered with some caution. All the forces are still widely 
based on phenomenology. Nevertheless the results indicate that very likely the effects will 
not be distributed rather smoothly over the whole phase space but the large effects will be 
concentrated in special regions. Our study might therefore help to guide future experimental 
searches. 

As an example of such guidance we performed calculations for the p(d,pp)n experiment 
currently under way at KVI Groningen. This experiment covers a good fraction of the full 
phase space. It will be very interesting to see how the upcoming data for cross section, 
vector and tensor analyzing powers will be related to our predictions. 

We have just begun to understand the role which 3NFs play in a nuclear Hamiltonian. 
Our present study, though not complete, clearly shows the significance and importance of 
the 3NF effects in the breakup process, especially at higher energies. A natural extension 
of this study is the four-nucleon scattering that might reveal even more sensitivity to the 
underlying forces. The first steps toward that goal have already been made |fT0"l-|72| . 

In order to fully understand the Nd scattering dynamics at higher energies it is necessary, 
in addition to the inclusion of 3NFs, to formulate the equations in a relativistically covariant 
way. Also exact inclusion of the pp Coulomb force in pd scattering at energies above the 
deuteron breakup threshold is still unsolved and leaves a very uncomfortable theoretical 
uncertainty in the analysis of pd data. The solution of these two problems: the relativistic 
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formulation of 3N scattering equations and the exact inclusion of the Coulomb interactions 
above the deuteron breakup threshold will remove the remaining uncertainties and will make 
Nd scattering an extremely precise tool for testing nuclear interactions. It will also allow 
to compute exact 3N wave functions that can be used to study electromagnetic processes in 
the 3N system for a wide range of energies. 



ACKNOWLEDGMENTS 

This work was supported by the Polish Committee for Scientific Research (Grant No. 
2P03B02818 and 5P03B12320), the Deutsche Forschungsgemeinschaft (J.G.), and the NSF 
(Grant No. PHY0070858). One of us (W.G.) would like to thank the Foundation for Polish 
Science for the financial support during his stay in Cracow. R.S. is a holder of scholarship 
from the Foundation for Polish Science and acknowledges financial support. The numerical 
calculations have been performed on the Cray T90 and T3E of the NIC in Jiilich, Germany. 



15 



REFERENCES 



W. Meier, W. Glockle, Phys. Lett. B138, 329 (1984); B. H. J. McKeller, W. Glockle, 
Nucl. Phys. A416, 435 (1984). 

Proceedings of the International Symposium on the Three-Body Force in the Three- 
Nucleon System, Washington D.C, April, 1986, Eds. B. L. Berman, B. F. Gibson, 
Lecture Notes in Physics 260 (1986). 

W. Glockle, H. Witala, D. Hiiber, H. Kamada and J. Golak, Phys. Rep. 274, 107 (1996). 
V. G. J. Stoks, R. A. M. Klomp, C. P. F. Terheggen, J. J. de Swart, Phys. Rev. C49, 
2950 (1994). 

R. B. Wiringa, V. G. J. Stoks, R. Schiavilla, Phys. Rev. C51, 38 (1995). 
R. Machleidt, F. Sammarruca, and Y. Song, Phys. Rev. C53, R1483 (1996). 
R. Machleidt, Phys. Rev. C63, 024001 (2001). 

S. Weinberg, Phys. Lett. B251, 288 (1990); Nucl. Phys. B363, 3 (1991). 
U. van Kolck, Phys. Rev. C49, 2932 (1994). 

C. Ordonez, L. Ray, U. van Kolck, Phys. Rev. C53, 2086 (1996). 

T.-S. Park, K. Kubodera, D.-P. Min, M. Rho, Nucl. Phys. A646, 83 (1999). 
N. Kaiser, R. Brockmann, W. Weise, Nucl. Phys. A625, 758 (1997). 
N. Kaiser, S. Gerstendorfer, W. Weise, Nucl. Phys. A637, 395 (1998). 
E. Epelbaum, W. Glockle, U.-G. MeiBner, Nucl. Phys. A637, 107 (1998). 
E. Epelbaum, W. Glockle, U.-G. MeiBner, Nucl. Phys. A671, 295 (2000). 
M. R. Robilotta, Phys. Rev. C63, 044004 (2001). 

E. Epelbaum, H. Kamada, A. Nogga, H. Witala, W. Glockle, U.-G. MeiBner, Phys. Rev. 
Lett. 86, 4787 (2001). 

D. R. Entem, R. Machleidt, Phys. Lett. B524, 93 (2002). 

E. Epelbaum et al, |nucl-th/020TUB4. 



E. Epelbaum et al, to appear in Proceedings of the conference on "Mesons and Light 
Nuclei", Prague, July 2001 (|nucl-th/0109065| ). 



D. R. Entem, R. Machleidt, H. Witala, |nucl-th/0111033i 

A. Nogga, H. Kamada, and W. Glockle, Phys. Rev. Lett. 85, 944 (2000). 

M. Viviani, Nucl. Phys. A631, 111c (1988). 

A. Nogga, H. Kamada, W. Glockle, B.R. Barrett, frucl-th/0112026 . 



J. Carlson, R. Schiavilla, Rev. of Mod. Phys. 70, 743 (1998). 

R. B. Wiringa, S. C. Pieper, J. Carlson, and V. R. Pandharipande, Phys. Rev. C62, 
014001 (2000). 

J. Fujita, H. Miyazawa, Prog. Theor. Phys. 17, 360 (1957). 

B. S. Pudliner, V. R. Pandharipande, J. Carlson, S. C. Pieper, and R. B. Wiringa, 

Phys. Rev. C56, 1720 (1997). 

S. A. Coon et al, Nucl. Phys. A317, 242 (1979). 

S. C. Pieper et al, Phys. Rev. C64, 014001 (2001). 

K. Sekiguchi et al, Phys. Rev. C65, 034003 (2002). 

H. O. Meyer, private communication. 

K. Bodek, private communication. 

K. Sagara, private communication. 

H. Sakai et a/.,Phys. Rev. Lett. 84, 5288 (2000). 

N. Sakamoto et al, Phys. Lett. B367, 60 (1996). 



16 



[37] W. Tornow, C. R. Howell, R. C. Byrd, R. S. Peroni, and R. L. Walter, Phys. Rev. Lett. 

49, 312 (1982). 
38] R. V. Cadman et al, Phys. Rev. Lett. 86, 967 (2001). 
39] R. Bieber et al, Phys. Rev. Lett. 84, 606 (2000). 
40] K. Ermisch et al, Phys. Rev. Lett. 86, 5862 (2001). 
41] P. Hempen et al, Phys. Rev. C57, 484 (1998). 
42] W. Griiebler et al, Nucl. Phys. A353, 31c (1981). 

43] L. Sydow et al, Nucl. Phys. A567, 55 (1994); L. Sydow et al, Few-Body Syst. 25, 133 
(1998). 

44] H. Witala et al, Phys. Rev. C63, 024007 (2001). 
45] St. Kistryn et al, Nucl. Phys. A689, 345c (2001). 

46] D. Huber, H. Kamada, H. Witala, and W. Glockle, Acta Physica Polonica B28, 1677 
(1997). 

47] A. Stadler, W. Glockle, P. U. Sauer, Phys. Rev. C44, 2319 (1991). 
48] G. O. Ohlsen, Rep. Prog. Phys. 35 (1972) 717. 

49] G. G. Ohlsen, R. E. Brown, F. D. Correll, and R. A. Hardekopf, Nucl. Instr. and 

Methods 179, 283 (1981). 
50] S. A. Coon, W. Glockle, Phys. Rev. C23, 1790 (1981). 
51] S. A. Coon, H. P. Han, Few-Body Syst. 30, 131 (2001). 

52] For a short review see M. R. Robilotta, Few-Body Syst. Suppl. 2, 35 (1987). 

53] S. A. Coon, M. T. Pefia, Phys. Rev. C48, 2559 (1993). 

54] S. A. Coon, M. T. Pefia, D. O. Riska, Phys. Rev. C52, 2925 (1995). 

55] B. D. Keister, R. B. Wiringa, Phys. Lett. B173, 5 (1986). 

56] T.-Y. Saito, J. Haidenbauer, |nucl-th/0003064| . 

57] H. Witala et al, Phys. Rev. C52, 1254 (1995). 

58] J. L. Friar, D. Huber, U. van Kolck, Phys. Rev. C59, 53 (1999). 

59] D. Huber, J. L. Friar, A. Nogga, H. Witala, U. van Kolck, Few-Body Syst. 30, 95 (2001). 

60] M. R. Robilotta, H. T. Coelho, Nucl. Phys. A460, 645 (1986). 

61] A. Kievsky, M. Viviani, and S. Rosati, Phys. Rev. C52, R15 (1995). 

62] J. Carlson, V. R. Pandharipande, R. B. Wiringa, Nucl. Phys. A401, 59 (1983). 

63] H. Witala, W. Glockle, J. Golak, D. Huber, H. Kamada, and A. Nogga, Phys. Lett. B 

447, 216 (1999). 
64] P. U. Sauer, Prog. Part. Nucl. Phys. 16, 35 (1986). 
65] S. Nemoto et al, Phys. Rev. C58, 2599 (1999). 

66] A. Nogga, D. Huber, H. Kamada, and W. Glockle, Phys. Lett. B 409, 19 (1997). 
67] W. P. Abfalterer, et al. Phys. Rev. Lett. 81, 57 (1998). 
68] H. Witala et al, Phys. Rev. C59, 3035 (1999). 
69] J. Kuros-Zolnierczuk, PhD thesis, Krakow, 2001, 

[ht t p : / / zlj - www . if . uj . edu . pi/ pr j / 1 hf ew b o d . ht m| 
70] M. Viviani, S. Rosati, and A. Kievsky, Phys. Rev. Lett. 81, 1580 (1998). 
71] A. Fonseca, Phys. Rev. Lett. 83, 4021 (1999). 

72] F. Ciesielski, J. Carbonell, and C. Gignoux, Phys. Lett. B447, 199 (1999). 

73] J. Seagrave et a/., Ann. Phys. (N.Y.) 74, 250 (1972). 

74] G. Pauletta and F. Brooks, Nucl. Phys. A255, 267 (1975). 

75] M. Holmberg, Nucl. Phys. A129, 327 (1969). 



17 



[76] H. Catron et al, Phys. Rev. 123, 218 (1961). 



18 



TABLES 



TABLE I. Comparison of the nd total breakup cross section a^ ot with the value crf^ n evaluated 



according to Eq.(ll) at four energies of the incoming nucleon. 



Eiab [MeV] 


°tt [nib] 


<? n [mb] 


13 


167.5 


159.2 


65 


91.7 


87.3 


135 


56.0 


53.5 


200 


48.2 


45.9 



TABLE II. The maximal changes AO = max[<50(Ei, E2, 61, 02, ^12)] due to 3NF effects for the 
eight studied observables. For each energy the first row, labeled all, gives AO without constraints, 
the second row, labeled cr-cut, gives AO under the condition a^ 2N ^ > 0.01 sr -I^ e y , the third row, 
labelled cr + E-cut, gives AO under the additional restriction that E\^ > 15(3) MeV for Ei ao = 65, 
135, 200 (13) MeV. The last column gives the total number of studied configurations (covering all 
values of AO). 



Elab 


cuts 


Act 


AA X 


AA Z 


AAy 






AAyy 


AA;£ 2 


N 


[MeV] 




[%] 


[%] 


[%] 


[%] 


[%} 


[%] 


[%] 


[%] 




13 


all 


10 








7 


10 


22 


32 


19 


6147 




cr-cut 


10 








7 


10 


22 


32 


19 


6008 




a + E-cut 


10 














18 


17 


15 


1532 


65 


all 


23 


110 


33 


182 


88 


188 


194 


211 


7893 




cr-cut 


23 


97 


33 


156 


88 


188 


194 


205 


7323 




a + E-cut 


23 


97 


22 


156 


88 


188 


126 


130 


2002 


135 


all 


98 


202 


99 


327 


289 


357 


281 


577 


8253 




cr-cut 


98 


93 


58 


269 


233 


351 


166 


195 


7370 




a + E-cut 


98 


76 


58 


269 


233 


351 


154 


195 


1028 


200 


all 


133 


251 


205 


418 


458 


580 


702 


677 


8597 




cr-cut 


115 


81 


60 


238 


368 


275 


443 


230 


6638 




a + E-cut 


115 


59 


60 


238 


368 


275 


443 


230 


689 
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TABLE III. The same as in Table II for five observables in the p(d,pp)n breakup reaction at 
130 MeV. 



cuts 


Act 






AAyy 


AA XZ 


N 


all 


23 


85 


197 


181 


172 


970 


a-cut 


23 


85 


197 


181 


172 


970 


a + .E-cut 


21 


80 


100 


181 


172 


955 



20 



1000 




lab 



30 40 

[MeV] 



50 



FIG. 3. Comparison of the CD Bonn potential predictions (solid line) with the experimental 



data (full dots) |73-(fq| for the total elastic cross section o"fo( S . 
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FIG. 5. Comparison of theoretical predictions with and without 3NFs to the experimental data 
(full dots) |57| for the nd total cross section. The solid and dashed lines are the CD Bonn and CD 
Bonn+TM 3NF predictions, respectively. The open squares are the results for AYT8+Urbana IX 
and the circles for CD Bonn+TM'. 
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FIG. 6. 3NF effects in the nd total cross section (solid line), the total elastic scattering cross 
section (dotted line), and the total breakup cross section (dashed line), expressed as a percentage 
change (see Eq. [10| ) with respect to the total cross section as obtained with the NN interactions 
alone: CD Bonn+TM (top), CD Bonn+TM' (middle), and AV18+Urbana IX (bottom). 
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FIG. 7. Two-dimensional projections of phase space regions for Act in the Nd breakup at 
135 MeV. In the first row (all) phase space points are plotted without any restrictions, in the 
second row (cr-cut), and the third row (a + E-cut) the remaining points are plotted after the 
corresponding cuts have been performed. The three columns show the 9\ — 62, #1 — <f>u and 
Ei — E2 planes, respectively. The different colours express the Act values in percentage. 
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FIG. 8. Two-dimensional projections of phase space regions for AA yy in the Nd breakup at 
135 MeV. For the description see caption of Fig. |?j. 
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FIG. 9. Cross section in [mb MeV _1 sr -2 ] and analyzing powers for selected breakup configu- 
rations at 13 MeV. The light shaded band (blue) contains the theoretical predictions based on CD 
Bonn, AV18, Nijm I, II and Nijm 93. The darker band (red) represents predictions when these NN 
forces are combined with the TM 3NF. The solid line is for AV18+Urbana IX and the dashed line 
for CD Bonn+TM'. 
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FIG. 10. Differential cross section in [mb MeV sr 2 ] in selected breakup configurations at 65 
(top), 135 (middle) and 200 MeV (bottom). For the description of bands and lines see caption of 
Fig. | 
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FIG. 11. Nucleon analyzing power A y in selected breakup configurations at 65 (top), 135 (mid- 
dle) and 200 MeV (bottom) . For the description of bands and lines see caption of Fig. |j| 
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FIG. 12. Deuteron vector analyzing power Ajf* in selected breakup configurations at 65 (top), 
135 (middle) and 200 MeV (bottom). For the description of bands and lines see caption of Fig. ||. 
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FIG. 13. Nucleon analyzing power A z in selected breakup configurations at 65 (top), 135 (mid- 
dle) and 200 MeV (bottom) . For the description of bands and lines see caption of Fig. |j| 
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FIG. 14. Nucleoli analyzing power A x in selected breakup configurations at 65 (top), 135 
(middle) and 200 MeV (bottom). For the description of bands and lines see caption of Fig. ||. 



34 




<3 0.0 



S [MeV] 



FIG. 15. Tensor analyzing power A yy in selected breakup configurations at 65 (top), 135 (mid- 
dle) and 200 MeV (bottom) . For the description of bands and lines see caption of Fig. |j| 
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FIG. 16. Tensor analyzing power A xx in selected breakup configurations at 65 (top), 135 (mid- 
dle) and 200 MeV (bottom) . For the description of bands and lines see caption of Fig. [9|. 
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FIG. 17. Tensor analyzing power A xz in selected breakup configurations at 65 (top), 135 (mid- 
dle) and 200 MeV (bottom) . For the description of bands and lines see caption of Fig. |j| 
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FIG. 18. Theoretical predictions for the cross section in [mb MeV _1 sr -2 ] of the dN breakup at 
130 MeV for selected configurations. For the description of bands and lines see caption of Fig. ||. 
The points along the 5-curve displayed in the figures correspond to the ones accessible in the 
Groningen experiment. 
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FIG. 19. Theoretical predictions for the deuteron vector analyzing power of the dN breakup 



at 130 MeV for selected configurations. For the description of bands and lines see caption of Fig. 18 



39 




120 140 160 180 200 220 

T 




S [MeV] 



FIG. 20. Theoretical predictions for the tensor analyzing power A xx of the dN breakup at 
130 MeV for selected configurations. For the description of bands and lines see caption of Fig. Il8|. 
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FIG. 21. Theoretical predictions for the tensor analyzing power A yy of the dN breakup at 
130 MeV for selected configurations. For the description of bands and lines see caption of Fig. [l8|. 
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FIG. 22. Theoretical predictions for the tensor analyzing power A xz of the dN breakup at 
130 MeV for selected configurations. For the description of bands and lines see caption of Fig. Il8|. 
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